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, We consider the adjacency matrix A of a large random graph and study fluctuations of 

fT) ■ the function f n (z,u) — — Y^l=i ey A>{— u Gkk{z)} with G(z) = [z — iA)^ 1 . We prove that 

the moments of fluctuations normalized by n~ 1//2 in the limit n — > oo satisfy the Wick 
r*j ' relations for the Gaussian random variables. This allows us to prove central limit theorem 

<*S ^ for Tr G(z) and then extend the result on the linear eigenvalue statistics Trip(A) of any 

function tp : R —> K which increases, together with its first two derivatives, at infinity not 
faster than an exponential. 
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1 Introduction 



■ Random graphs appear in different branches of mathematics and physics (see monographs 

[H [12] and references there in). It is well known that they are closely connected with the 
theory of random matrices, since there is one to one map between graphs with n vertices 
I/"") ■ and their adjacency matrices (recall that by the definition the entries a^j of the adjacency 

matrix are 1 if the vertices i and j are connected and a%j = otherwise). Commonly, the 
set of n eigenvalues of the adjacency matrix is referred to as the spectrum of the graph. The 
limit when the dimension of the matrix n (the number of the vertexes of the graphs) tends 
to infinity provides a natural approximation for the spectral properties of random graphs. 

One of the classes of the prime reference in the theory of random graphs is the binomial 
random graph originating by P. Erdos (see, e.g. [12]). Given a number p n £ (0, 1), this family 
1 of graphs G(n,p n ) is defined by taking the set of all graphs on n vertices as the space of 

events with probability 

P(G)=p e n ^(l-p n p-< G \ (1.1) 

where e{G) is the number of edges of G. Most of the random graphs studies are devoted to 
the cases where p n — > as n — > oo. 

Ensemble of random symmetric n x n adjacency matrices A corresponding to (jl.ip can 
be represented as A = {ajj}" J=1 with an = 0, and i.i.d. 

1, with probability p n , . . 

0, with probability 1 — p n , 

For any measurable function / we denote E{/(j4)} the averaging with respect to all random 
variables {p4j}x<i<i<n an d 

Var{/(^)} := B{\f(A) - E{f(A)}\ 2 }. (1.3) 
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The normalized eigenvalue counting measure of A is denned by the formula 

N n (\)=n-H{j:\ ( f ) <\}. 

The ensemble of adjacency matrices (ll.2p is a particular case of the random matrix theory, 
where the limiting transition n — > oo is intensively studied during half of century since 
the pioneering works by E. Wigner [23]. Spectral properties of random adjacency matrix 
(ll.2p were examined in the limit n —* oo both in numerical and theoretical physics studies 
[3 O [18l UHl [20]. There are two major asymptotic regimes: p n 3> 1/n and p n = 0(l/n) 
and corresponding models can be called dilute random matrices and sparse random matrices, 
respectively. The first studies of spectral properties of sparse and dilute random matrices 
in the physical literature are related with the works [19] . [20] . [18] . where equations for the 
limiting density of states of sparse random matrices were derived. In papers [18] and [10] a 
number of important results on the universality of the correlation functions and the Anderson 
localization transition were obtained. Unfortunately these results were obtained with the non 
rigorous replica and super symmetry methods. 

On mathematical level of rigor the eigenvalue distribution of dilute random matrices was 
studied in p3j. It was shown that the normalized eigenvalue counting measure of (np n )~ l l 2 A 
converges in the limit np n — > oo to the distribution of explicit form known as the semicircle, 
or Wigner law [23]. In the paper [5] the adjacency matrix of random graphs (jl.ip with 
p n = pn~ l was studied. It was shown that for any m there exist non random limiting 
moments lirrin^oo n Tr A™ and these moments can be found from the system of certain 
recurrent relations. The results of [5] was generalized to the case of weighted random graphs 
in [15] , where the resolvent of the adjacency matrix was studied and equations for the Stieltjes 
transform g(z) of the limiting eigenvalue distribution were derived rigorously (note, that the 
same equation for gaussian weights were obtained in [19], [20], [18] by using the replica and 
the super symmetry approaches.) It was shown in [15] that to prove the existence of the 
limit lirrin^oo g n (z) = g(z), where g n (z) is the Stieltjes transform of the normalized counting 
function N n (X) 

f dNJX) 

we need to study the behavior of the function 

1 - 

f n (z,u) = -J2^' uGkk(z \ (1-5) 
k=l 

where 

Gjk(z) = (z-iA)j£, 

The function f n (z,u) is defined for any u,z such that "Stz 
important for us that 

n 

\\G\\< \dtz\-\ Yl \°ij\ 2 = (GG*)u < \\G\\ 2 < |3te|- 2 (1.7) 

3=1 

K(Ge,e)Kz>0, Ve G M n ^ \e~ u( - Ge ' e) \ < 1, if u^z > 0. 

Here and everywhere below ||^4|| means the operator norm of the matrix A. 

The following theorem (proven in [15]) gives us the limiting properties of f n (z,u) of (|1.5p 



(1.6) 

7^ 0. In what follows it will be 
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Theorem 1 Consider the adjacency matrix with p n = p/n. Then for any u,z such 

that vMtz > we have: 

(i) the variance of the function f n (z,u) defined by ( li.51) vanishes in the limit n —* oo: 

\sx{f n (z,u)} <C/(Rz) 2 n, (1.8) 

(ii) there exists the limit 

lim E{f n (z, it)} = f(z, u), \B{f n (z, it)} - f(z, u)\ < Cu 1 / 2 /\Rz\n l l 2 (1.9) 
n— >oo 

(m) tf we consider a class H of functions which are analytic in z : $lz > and for any fixed 
z : > possessing the norm 

||/(z)||=max^dl, (1.10) 
i/ien i/ie limiting function is the unique solution in TC of the functional equation 

f(z, u) = 1 - uV2 e -P /"°° d / l(2 )^ exp{-^ + p/(s, u)}, (1.11) 

JO y/V ■ 

where is the Bessel function 



One can easily see that 

1 " 1 

- VE{G fcfc (z)} = — E{Tr G(z)} = B{ig n (-iz)}, 



u-0 k=1 



where g n (z) is the Stieltjes transform (jl.4p of the normalized counting measure N n (X). Hence, 
Theorem [T] implies that for any z : Qz ^ 

lim E{\g n (z) - E{g n (z)}\ 2 } = 0, (1.13) 

n— >oo 

i.e., the fluctuations of g n {z) vanish in the limit n — > oo. And (11. 9j) implies that 

(1.14) 



= hm E{g n (z)} = -—f(z,u) 

ti—kx> Oil 



u=0 



Since the Stieltjes transform uniquely determines the measure, it follows from Theorem [T] 
that there exists the weak limit N(X) of the normalized counting measure N n (X) and the 
Stieltjes transform g{—iz) can be obtained as the first derivative of the solution of (|l.lip . 
Using Theorem Q] it is not difficult to obtain the asymptotic expansions for g(z) with respect 
to z~ k . Since it is well known that the coefficients of this expansion are the moments of the 
limiting normalized counting measure of eigenvalues, we obtain the recurrent formulas for 
the moments. Besides, constructing the asymptotic expansion of g(z) with respect to p k , it 
is easy to show that this expansion is convergent for p < 1. Since in the case a™ = 0, 1 the 
coefficients of this expansion are rational functions of z, we can conclude that the limiting 
spectrum is pure point and consists of the spectra of finite blocks only. 
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Results of [15] described above can be viewed as the analogs of the Law of Large Numbers 
for linear eigenvalue statistics 



jV n [(p] = J2<p(\i) = Tr<p(A) (1.15) 
i=l 

corresponding to continuous test functions. Indeed, it follows from (jl,13p - (|1.14p that for 
any continuous test function there exists 

lim n^NrM = I V>WdN(\), 

n—>oa J 

where N is the limiting normalized counting measure of eigenvalues. In the present paper we 
consider the central limit theorem, the second element of the standard probabilistic analysis 
of linear statistics. Similar questions for other ensembles of random matrices were studied in 
[31 [HI [HI [HI [211 [22] . Note, however, that for almost all ensembles studied in the random 
matrix theory, like the Wigner ensemble, the Marchenko-Pastur ensemble, the matrix models, 
etc the variance of linear statistics for smooth functions is bounded (see [21 [31 [TT], [131 Q51 ED 
22J). Thus, for these ensembles, one expects the Central Limit Theorem to be valid for 
statistics themselves, i.e., without an n-dependent normalization factor in front. This has 
to be compared with the case of i.i.d. random variables with finite second moment, where 
the variance of linear statistics is always of the order 0(n), n — > oo and the Central Limit 
Theorem is valid for linear statistics divided by n 1//2 . As we will see below this is the case 
also for the ensemble of sparse adjacency matrices (11. 2\i with p n = p/n. 

The aim of the present paper is to study the fluctuations of linear eigenvalue statistics for 
different classes of test functions. Following the method of [15] we study first the functions 
f n (z,u) (defined in (jl.5p ) and prove that its fluctuations converges in distribution to the 
complex Gaussian random variables. 

Define the m-th generalized moment of the fluctuations of f n (z, u): 

{m / n _ u .q r z .\ 
n e ? 
3=1 \k=l 



n 



m/2 E 



Hfn(Zj,Uj)\, MZi^O. (1.16) 
0=1 



Here and below for any random variable £ we denote 

I = e - eu} 

Theorem 2 Consider the adjacency matrix \1.2\i withp n = p/n. Let M m n (zx,ux; . . . ; z m , u m ) 
(m = 2, 3, . . . ) of §1.16}) be the "moments" of the fluctuations of f n (z, u) of 1 1.5(1 . Then for 
any m > 2 and z\ , . . . , z m : $tzj > there exists 

M m (zx,ux;...;z m ,u m ):= lim M mn (zx,ux;...;z m ,u m ). (1.17) 

n— >oo 

Moreover, the following recursion equations hold: 
M m { 

Zx j ^1 1 ■ • • > Zm j U m ) 

rn 

= } y M 2 (zx , ux ; Zj , nj)M m -2 (z 2 , u% ; . . . ; Zj-i, u 3 -x; z j+ x, %+i; • • •;z m ,u m ). (1.18) 

3=2 
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Theorem [2] can be used to prove the central limit theorem for fluctuations of the trace of 
G(z) of (jl.6p . Indeed, if we denote 

M* tn ( Zl , ...,z m ):= n~ m / 2 B jlr G{z x ) . . . TV G(z m ) j , (1.19) 

then it is easy to see that 



Qm 

M*(zi, ...,z m ) = - 5 — M m (zi,ui; . . . ;z m ,u m ) 

OU\ . . . ou m 



ui=--=u m =0 



Since M m (zi,ui; . . . ; z m , u m ) are evidently analytic in each U{ in some neighborhood of Ui = 
and bounded uniformly in n for any fixed Z\, . . . , z m (3?Zj ^ 0) (see Lemma [1] below), we pass 
to the limit n — > oo in the above relations and obtain the following theorem: 

Theorem 3 Let G(z) be the resolvent S1.6\) of the sparse adjacency matrix il.tfy with p n = 
p/n. Then for any m > 2 and z±, . . . ,z m : Wtzj > there exists 

M* m (z 1 ,...,z m ):= lim M* {z\; . . . ; z m ) (1.20) 

n— >oo ' 

and the following recursions hold: 

m 

M^(zi, ...,z m ) = ^M 2 * (zi,Zj)M^_ 2 (z 2 , . . .,Zj- 1 ,z j+1 , ...,z m ) (1.21) 

i=2 

o 

Theorem [3] by a standard way implies the central limit theorem for v n (z) = n-^Tr G(z). 
Indeed, if we put in (jl.20p - (jl.2ip z\ = z 2 = ■■• = z m = z, then Theorem [3] yields that there 
exist limits of all moments of the complex random variable v n (z) and 

M 2m (z) := hmEUn-^Tr G(z)) 2m } = (2m - l)\\{M 2 {z)) m 

This means that v n (z) converges in distribution to a complex Gaussian random variable with 
zero mean and variance M 2 (z). 

It is possible also to derive from Theorem [3] the central limit theorem for the linear 
eigenvalue statistics of any function ip which grows not faster than an exponent at infinity 
and possesses two derivatives with the same property, i.e. there exists a constant c > such 
that ip, ip', ip" G L 2 (M., cosh~ 2 (cA)). Here and below 

L 2 (R,w(X)) = (/ : f \f(X)\ 2 w(X)d\ < oo) (1.22) 



Theorem 4 Consider the adjacency matrix hi. 2^ with p n = p/n and take any function (p 
which possesses two derivatives such that ip, ip', ip" G L 2 (M, cosh~ 2 (cA)) with some constant 

o 

c > 0. Then the random variable n~ l ' 2 J\f n [<p\ converges in distribution to a Gaussian random 
variable with zero mean and variance VUp] := lim Var{n _1 / 2 A/" n [i^]}. 

n^oo 

It is clear from the above discussion that Theorem [2] plays a key role in the paper, because 
Theorems [3] and [4] are in fact corollaries of Theorem [2j The proof of Theorem [2] is based on a 
version of the cavity method which has been used many times for proving different limiting 
relations of statistical mechanics and random matrices. The idea is to compare the behavior 
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of the object function (e.g. free-energy, resolvent, etc.) for the complete system of random 
variables of the problem and the one with some subset of random variables replaced by 0. 

Let us try to explain the connections among the lemmas and propositions which are 
necessary for the proof of Theorem [2j The proof should be seen as a logical sequence of the 
following steps: 

• We prove first bounds on M mi( ,(zi, U\\ . . . ; z m , u m ) uniform in (z\, u\\ . . . ; z mi u m 

C > 0) (see Lemma [TJ. One uses the norm estimates of the martingale theory (Propo- 
sition [1]) , identities for the resolvent and the cavity method consisting in studying the 
difference of the resolvent of the full matrix and the same matrix without the first line 
and the first column. 

• To prove the convergence of the variance of the sums of exponentials we need to gener- 
alize Theorem [T] and to show the existence of the limits of exponentials multiplied by 
some entire functions (cf Lemma [2] and Theorem [I]) . The proof of Lemma [2] is based 
on the relations for some functions of A given by Proposition [2j 

• Lemma [3] proves the self averaging properties and the existence of the limits for the 
terms which will appear in the proof of CLT. 

• Finally we prove that the " moments" (11 . 16f) as functions of it, satisfy the linear integral 
equations with the kernel defined in terms of the function / of (see (|2.54p ). Since 
we are able to prove that these equations are uniquely solvable for Jtz > Mq with some 
fixed Mq, we finish the proof of Theorem [2j 



2 Proofs 

We start from the lemma which gives bounds for M mn . 

Lemma 1 For any m £ N and z\ , . . . , z m : ^tzj > there exists a constant C m such that 
uniformly in U\ . . . , u m > 

\M mtn (zi,u 1 ;...;z m ,u m )\<C m (2.1) 

The proof is based on the martingale property of the sequence of averages of the functions 
of the random matrix A with respect to its rows or columns. The sequence is ordered with 
respect to the index of the rows and the proposition below is based on the sequence of the 
conditional expectations like in the proof of self-averaging of the free-energy for disordered 
systems. 

Proposition 1 Let £ a , a = be independent random variables, assuming values in 

R ma and having probability laws P a , a = l,...,v and let <3? : M mi x • • • x W 1 " — > C be a 
Borelian function. Set 

= J Za,Za+l,---,Zv)P a +l(dZ a+ i)...P v { ( % l ,) (2.2) 

so that <$> u = <3?o = E{<1>}, where E{...} denotes the expectation with respect to the 
product measure P\ . . . P u . 

Then for any positive p > 1 there exists C' p , independent of v and such that 

V 

E{|<D - E{<D}| 2 f} < Cy- 1 V{\$a ~ $a-i| 2p }. (2.3) 
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Moreover, if for every a = 1, . . . , v there exists a ^-independent : R mi x • • • x 
such that 



then 



E{|$-^ a )| 2p }<G<oo, a = l, 
E {|$-E{$}| 2 P} < 2C' p Cv p . 



-> C 
(2.4) 

(2.5) 



Proof. The proof of (|2.3p is given in [6]. Hence, we show only how to derive (|2.4p from 
(1231) . It follows from and that the integrals of with respect to P a+ i...P w and 
P a P a +i...P u coincide and we obtain 

E{|*« - * a _ipP} < 2 2p - 1 (E{|($-^( a )) a _ 1 | 2p } + E{|($-^( Q )) a | 2 P}) 
< 2 2p E{|$ - ^ (a) | 2p }. 

This and ([23]) prove (j23]h □ 

Proof of Lemma [7] The Holder inequality yields 



\M m ^ n (z 1 



<n m / 2 J]E| 

7 = 1 ^ 



m\ l/m 



Hence, it suffices to prove the bound for the r.h.s. of the above inequality. For this we use 
Proposition Q] for the function $ = nf n (z,u) with f n {z,u) of (|1.5p . According to (|2.5p and 
the approach of the cavity method for our purposes it is enough to choose the functions 
independent from £j = := {an, • • • , a>u-i, 0, Qii+i, . . . , aj n ) and prove (|2.4p . Set 



.4 



, G®(z) = (z-iA®)-\ 



n 



ciij =0J=l,...,n 

/W(z,u) :=J> 



(2.6) 
(2.7) 



By the symmetry reason it suffices to prove (|2,4p for i = 1. We use the representations: 



G«(z) = CS?(z) 



^GWaW), 



z+ ( G (i) a (i),a(i)) 



(2.8) 



where = (0, ai2, . . . , a\ n ). The inequality \e x — e y \ < \x — y\max{\e x \,\e y \} and (|1.7j) 
imply 



fc=2 



fc=2 



(2.9) 



(G( 1 )a«) jfc (G( 1 )a( 1 )) fc (GW^aW,^ 1 )^)^ 1 )) 



| 2 + ( G (i) a (i) 5a (i))| 



Iz+CGWoW.aW)! 



But the spectral theorem yields 



{G^{z)a^,G^{z)aW) = J2 



|(^),aW)p 
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where A^i/j^ = A^ (j) - Thus, since by (H2J) ^z^(G^a^, a^) > 0, we have 



(2.10) 



Inequality (|2.4|) for our choice of $ and $W follows from (|2.9p and f|2.10j> . □ 

In the proof of Theorem [5] we will replace sometimes M mn by the moments independent 
of {aij-}y =2 . Set 



AfW (zi,Mi;...;z m ,ii m ) := n m/2 E 



A o(1) 



(2.11) 



with fn °f (|2.7p . Note that (|2,9p yields that for any m £ N and z\ , . . . , z m : ?R.Zj > there 
exists constants C m , C' m such that uniformly in u± . . . , u m 

\M m ^ n (zi,ui; . . .;z m ,u m ) - M^ n (z 1 ,ur, . . . ;z m ,u m )\ < C m n~ 1/2 , 
\M$ n {z 1 ,u 1 ;...-z m ,u m )\ <C' m . (2.12) 

To study the behavior of some functions, depending on {aij}"_ 2 , we use the proposition: 

Proposition 2 Let Ei{. . . } be the averaging with respect to {aifc}fc = 2- Then we have for 
any u, v > and "tftz > 2 



(2.13) 



Y^GVauav + otv^GVauaifX E\ /2 {\r v \ 2 } < Cvn^ . 



Moreover, denoting Z = z + (G"0",o^), we /iaue 
f n i ( n 

^ i=2 J ^ i,jfc=: 



+0 



e ""33 (e v "j fe 

j,k=2 

(u + u 2 )e u ' 



l)oi * ^ (2.14) 



Proof. Note that since |9ffcz| 1 < 1/2, everywhere below we will replace 1 by a constant. 
We need below the trivial bounds: 

\e a -e b \ < |o-6|max{|e a |,|e fe |}, \e a - e b - (a - b)\ < \a - b\ 2 max{|e a |, \e b \}. (2.15) 

The first bound and the second line of (jl.7p combined with (|2.8p imply 

(i) 



| e -KGWaW,a«)_ e -„E fe G«a lfe | < y 
Averaging the square of the bound we obtain 



EM r 
G fc 1 fc 2 a ifci a ifc2 



(2.16) 



Ci^ 2 ^ (i) 2 



(1) 1 2 , C 2 V 2 



(1) ^(1) c 3 « 2 



n- 



kifa 



n° 



ST MV Mi) , 



ki,k2,k-3 



n' 



k\,k2 



< 



C A v 2 



n 
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Here we used the bounds valid for any matrix A: 



E^ 



< n 



Ei^*i^ nl/a (Ei-^*i a ) ^ 1/2 ihi- ( 2 - 17 ) 



To prove (|2.14p we show first that 



^^{(expl-^} - exp{-uGg> + £(G$) 2 a lfe /Z} 
i=2 



i=2 



< Cn~ 1 ' 2 e u . (2.18) 
The second inequality of (12. 15ft and the bounds that \e aj \ < 1 and |e fcj | < e"/' 3 * 2 ' 3 < e u yield 



n n n 

£ Bl {^-e*} < + ^1 {|(^ - J>) ~ (a* - fy)|} 



3=2 



3=2 



i=2 



< 



J>ite-M + e«J> 1 {|a,-6 i | 2 } 

3=2 j=2 



Then, similarly to (|2.16p we have 



E E i{ a i-M 

3=2 



E E n E ^ik^jfa 01 *! 01 ^^ f 



i=2 v fci^fc. 



< uE 



1/2 

i 



E (G (1) G^) klM a lkl a lk2 \b{ /2 {Z- 2 } < Curi' 1 / 2 . 



Here we used also that in view of (jl.7p $tZ > 3ftz > 2. Moreover, similarly to (|2.16p . we 
obtain 

E a {\ aj - b 3 \ 2 } < uW Yl GflGfla lkl a lk2 Gf k \Gfla lk3 a lki } < Cu 2 n 2 . 
Summing with respect to j, we get f|2. 18j) . Besides, we have 

n 

E (ex P {-uG« + ^(G«) 2 W^} " ex P {-uG«} 

3=2 

n co n ,mf\~^ //^i( 1 )\2„ \m 



E 

j=2 m=l 



uG (i) ^ u (G - fc j a lk 



E e ~" E 



m\Z r 



m=l 



m\Z r 



where the remainder term r n admits the bound 

n 

E i{wi < E E Ei {i g £i 2 i g Si Wi fe } E 



u m (E k \G 



(l)|2\m-2 n .,2„u 



j=2 fci^fe 



m=2 



T jfc I ^ , Cu 2 e h 
2(m-2)\Z r - 



< 



n 



The averaging here is similar to (|2.16l) . Thus, we have proved (|2.14p D 
Set (cf ([TT5]I ) 



c 



fc+1 



fc=0 



fc!(Jfc + l)! 



-2zC 1 / 2 Ji(2< 1 /2 ) 



(2.19) 



Below we will need the following properties of 



sup 1^(01 < Ji(r), \J X {Q\ < ICKI + Jxdd)), sup 1^(01 < {l + Jx{r)) (2.20) 
|C|<r K\<r 

The following lemma is the analog of Theorem [T] for the function which will appear in the 
proof of Theorem [3j 

Lemma 2 For any u > 0, v £ C ; $tz > 2, and J\ of U.1B the random variable 



Vj, n (z,u,v) = n" 1 ^ (T uGkk Jx(vG%) (2.21) 
j,k=l 

possesses the property: 

Var{V J>n (z, u, v)} < n^qiu, \v\)(l + J?(\v\)) (2.22) 
where q(u, v) is a fixed polynomial. Moreover, there exists 

Vj{z,u,v) := lim E{Vj, n (z,u,v)}. (2.23) 

n— >oo 

and 

\r Jin (z,u,v)\ := | V Jtn (z,u, v) -Vj(z,u,v)\ < CrT x l 2 {l + Ji{\v\)). (2.24) 
Proof. According to Proposition Q] to prove (12.221) it is enough to prove that 



<q 1 (u,v)(l + J 1 (\v\)) (2.25) 



with polynomial q±. Then q = q\. In view of the second bound of (|2.20p . (jl.7p . (|2.8p . (|2.9p . 
and (|2.10p we have 



A (i,i) 



E 



< Cu\v\(l + Ji(H)) ^ - 



< CuM^M) 1 , \> < Cu\v\(l + Ji(\v\)). 



| z + ( G (i) a (i) ja (i))| 
Moreover, by the third bound of (|2.20p . we can write 



A (l,2) . = 



e- uG ^ (MvGlA-JMGl) 



(i) 



(1)^2^ 



< 



E' 

+C| W | 2 (1 + JM)) £ - G &l(IG fcJ l 3 + |Gg 



,(1)2 



(1) |3^ 



(2.26) 
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Then denoting Si the first sum in the r.h.s., we have in view of the first line of (12. 8|) . (II. 7|) . 
and (HTTP]) .: 



4f |z+(GWa(i),a(i))| v GjfcU) + ^' fc W 



To estimate £2 ~~ the second sum in the r.h.s. of (|2.26p we use that for any matrix M if we 
consider the matrix = {|M|f ,-}^- =1 , then 



I \M® 1 1 < sup ( £ \ Mij | 2 ) V2 sup ( ^ |Mi 



1/2 

|2 ^ <||M|| 2 . 



(2.27) 



Hence, the matrix with entries |Gfcj| 2 has the norm bounded by ||G|| 2 < |3?z| 2 . Then (12.2TH 
and (|2.10p imply for S2: 

So < i^i-iy K^WaW^IK^WaW^l f, c , ,2 , |C (D L , |2 
S2 " 2^ |^ + (G(i)a(i),a(i))| 



< 2|^|" 



| z+ ( G (i) a (i)yi))| 



< c. 



Thus, we have proved (|2.25p and so (|2.22p . 

To prove (|2.23p - (|2.24p it suffices to prove that for any m > 2 there exists 

V m (u,z) = lim eU^Y e- uGkk G^X = lim e( V e'^GfX, 



and 



E {E e_MGllG ^}- y - 



(u,z) 



< Cm(l + tt)/n 1/2 



(2.28) 



To average with respect to a*- 1 ) we use the second and the third line of (|2.8[) and the formulas: 



x ^m— 1 



(m-1)! 



l-u 1/2 / dv 



Ji(2yfwv) 



exp{— R v}, 



(2.29) 
(2.30) 



which are valid for any ?R.R > and u € C. Then we get 



T m {u) := E £e-«^G^. =EEi / d VlJ -± 



00 „,m— 1 

«1 -v x aW.aM)) 



t^EEi 




dv\dv2- 



JO 



x (m-l)! 
<- 1 l 7 1 (2^i) „ 
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* c u 2 ( m ~ I)' 



K+^ 2 )(2+(G( 1 )a( 1 ),a( 1 ))) 



EEi ^(G«a«) 



(m — 1)! 



00 /-oo 



JO 



-n^EEi ^(G«a«)5 



1 v^Jxii^n) 

dv\ dv2 ;=— ; — — e 



(2.31) 

-(v 1 +v 2 )(z+(G < - 1 '>a < - 1 \aW)) 



fv 2 {m ~ I)' 
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Using (I2.13|) and averaging with respect to {aii}, we have 

h,m = / n d Vl rm 1 _ n! e-^EEi|exp{-z;i^4 1) ai,i 

dvi-^——e~ VlZ ^T\ (l - - + -e~ VlG '' } ) + 0(n _1 ) (2.32) 
im — ill J-J-V n n J 



oo 





i 


(m — 


1)! 




•l 


v?~ 




(m — 


1)! 




•l 


v™~ 




(m — 


1)! 




•l 


v™~ 




(m — 


1)! 



roc u m ~ l 

Jo {m-l)l 



where 



\n,m\ < Cmn 1/2 , 

and we used first (|2.9j) — (|2. 10|) to replace fn\z,v\) by f n (z,v\), and then (jl.9p to replace 
f n (z,vi) by f{z,vt). Similarly 



'2,m 




/0 JO 

Moreover, using (|2.14p and (|2.13p . we obtain 



Jvi(m - 1)! 





-i 


(m — 


1)! 




-i 


<~ 




(m — 


1)! 




-i 


v™~ 




(m — 


1)! 



i.m 



/ ! " - / av ^TZ 7TF ^ U jk> "I* ex P\ "1 2^ ° W f ^ r 3,m 

j,k I 

| £ e -* G $ (G^)™ exp{-p + 

/' / ^i /' 1 ^ e-^ z e- p+ ^ z ^T m (vi) + r 3 , m , (2.34) 
(m — 1J! 

\r 3 , m (u)\ < Cmn- 1 / 2 . 
Here we used also the relation 

E /I Ee -..<' (G a. r } = r mW + o(2). 

L 3,* J 

which can be proved, similarly to (j2.25p . Repeating the argument used, for I^^rni 

we obtain 



<4,rra 



o Jo 



U2(m — 1)! 



k 4 , m (u)| 2 < Cmn -1/2 . (2.35) 

Collecting the above relations, we get in view of (|2.31j) the equation 

T m (u) = <f m (u) + K m (T m )(u) + r m (u), 
\r m (u)\ < Cmtl + y^r 1 / 2 , 
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where the function ip m (u) is defined by the r.h.s. of (I2.32|) and (|2.33)) and the integral operator 
K m is defined by the r.h.s. of (|2.34[) and (|2.35j) . It is easy to see that for "Stz > 2 the operator 
norm in the Banach space of the functions with the norm (ll.lOp satisfies the inequality 

\\K m \\ < q < 1. 

Hence, we get (|2.28p . Then summing with respect to m and taking into account the bounds 
for the remainder terms, we obtain (|2.23p . □ 

The next lemma is a technical one. We will use it in the proof of Theorem [2] below. 



Lemma 3 Set 

n 

DW{z,u) := n(f n (z,u) - f£\z,u)) = e~ uG ^ + £ (V mG « - e~ uG ^ 

1=2 

Then for $tz > 2 we have 

Var |Ei | e -'" G "( z )JJ < n -\ Var {Ei {d«(z,u)}} < e^gi^rT 1 , 

Var |Ei U- ulGll ^D^(z 2 ,u 2 )\\ < e U2 q 2 {u 1 ,u 2 )n- 1 . (2.37) 

with polynomial qi,q 2 . Moreover, if we denote 

V n ( Zl , Ul ;z 2 ,u 2 ) = Covx { e -«iGii(*0 D (i) ^ ^ j f (238) 

where Cov 1 {F 1 ,F 2 } := Ei{FiF 2 } - Ei{Fi}Ei{F 2 }, then there exists 

V(zi,ui;z 2 ,u 2 ) = lim V n (zi,ui; z 2 ,u 2 ) (2.39) 

n— >oo 

and for any fixed Zi,ui; z 2 ,u 2 

\V(z 1 ,u 1 ;z 2 ,u 2 ) - V n (z 1 ,ui;z 2 ,u 2 )\ < q 3 (ux, u 2 )e" 2 n~ 1/2 , (2.40) 
with polynomial qz- 

Proof of Lemma [3J The first bound of (|2.37p can be proved similarly to (|2.3ip - (|2.34j) . 
Indeed, according to (12.30P we have 

Zb(u) := E 1 L~» G A = l- M i/a El | jf W efo Jl( ^ e^+C^)»W^)))l, 
Then, averaging with respect to similarly to (|2.32p . we get 

T (U) = 1-^1 dv Jl ( 2 ^\ -vz e -V+PJM + r(); E{ | ro |2 } < C/n . 

Jo 

To prove the second bound of (|2.37p we use (|2.14p . which gives us that Ei{L)W - e-" G n } 
coincides with the r.h.s. of (|2.14p . Then (|2.30p for u = i(G^) 2 u applied to the r.h.s. of 
([2441 yields: 

n Jo ^ 31 Vv 



3,1 



+0{q(u)e u /n 1 ' 2 ) = |/ cfc w -1 ^^, u, u^e"^^^ + O^/n 1 / 2 ) 



o 



13 



with Vj n of (|2.2ip . Now the second inequality of (12.370 follows from Lemma if we use 
(|2.30p to integrate the bound for rj tn (z,u,uv) of (|2.24[) with respect to v. The third bound 
of (I2.37P follows from the first and the second one. 

Relations (|2.39p - (|2.40p can be proved if we repeat the argument (|2.3ip - (|2.35p and then 
apply Lemma On 

Now we are ready to prove Theorem [2j 

Proof of Theorem [1 Fix z%, . . . , z m such that > 2, i = 1, . . . m. We find first M 2 ^ n . 
Using the symmetry of the problem and Lemma [T] it is easy to see that 



.1/ 



0-UlGu(2 1 )°( 1 ) 



2,)i 



o-«iGii(«i) »W, . 

D {Z 2 ,U 2 ) 



nEjg " f n (z 2 ,u 2 ) | +E|e 

= T 1 + V n (z 1 ,ui;z 2 ,u 2 ), 
where V n {z\, u\\ z 2 , u 2 ) is defined in Lemma Relations (|2.8p . (|2.30p . and (|2.13p yield 



(2.41) 



Tl 



-nu 



1/2 



dv 



-nu 



1/2 



00 , J x {2y/^v 
dv- 



o(l) 

^"E /„ (z 2 ,u 2 )e 



o(l) 



(2.42) 



^P^E j / n ' (z 2 , n 2 ) n e -^«« + r . 



with r„ of (|2T3|) 

o(l) 

Since /„ (22^2) does not depend on {aij}™ =2 we can average with respect to ' and 
similarly to (|2.16p obtain 

|Ei{r„}| < C(v + v 2 )/n. 

We used that (|2.17p and the first bound of (jl.7p for The bound, the Schwarz inequal- 

ity, and Lemma Q] yield 



o(l) 



<n- L \>Rz\- 2 VUf n (z 2 ,u 2 )\ \<C(v + v 2 )/n 



(i) 

E< i f n i z 2,u 2 )r l 



Then, integrating with respect to v (recall that < 1) and averaging \\ k e 
over {aik}, we get similarly to f|2.32|) : 



Tl 



-nu 



1/2 



0(1) 



dv -^ UlU ' e-^E j/ n ' (z 2 , u 2 ) (ft e~<^J j + 0({u/n)^) 



-nu 



1/2 



Writing f£\zi,v) = E{f n x> (z x ,v)} + f n (zt,v), we have 



r(l) 



0(1) 



-nu 



1/2 



dv 



^H#!) e -..-, E (/!; , (22 , 112)ef 



(*!,«)}+*>/„ (zi,v) I , (2) 



1 r ^ 

I ' Ti 



O(l) 



O(l) 



n / du^fuijUj^E^/n (z 2 ,u 2 ) [pfn (z 1 ,v) + 0[(f n (z 1 ,v)) >+rX 



0(1) 



dvK n (ui,v; zx)M 2 , n (zi,v; z 2 ,u 2 ) + r£ 



(3) 



(2.43) 
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where 

K n (u 1 ,v;z 1 ) :=- ^^M e -^'H} 

V v 

r°° fo(i) / 0(1) 

r (3) = r ^+nj dvK n (u 1 ,v;z 1 )Blf n {z 2 ,u 2 )0\{f n ( Zl ,v)f 

< Cn-^u 1 / 2 . (2.44) 

The last bound follows from (|2.12p . 
Thus, we obtain that 

M 2 ,n(zi,ui;z2,u 2 ) = / dvK n (ui,v;z 1 )M^l(z 1 ,v,Z2,u 2 ) + 

Jo 

V n (zx,ur,z 2 ,u 2 ) +r£\z 1 ,u 1 ;z 2 ,u 2 ) + 0(n" 1/2 ), (2.45) 

where V n (zi,ui] z 2 , u 2 ) is defined in (|2.38[) . Besides, using (|1.5p and the inequality |c7i(x)| < 1, 
we obtain that uniformly in u, v > 

lim K n (u, v; z) = _^i/2 ^i( 2 V^) e ~zv- P exp {_ E {p/(^ =: K ^ v . z ) sr z > 2 , 

and 

\K n (u, v; z) - K{u, v;z)\< Cuv^^e-^rT^I 2 . 

Using the above bounds to replace K n by K in (|2.45p and (|2.40p to replace V n by V, we can 
write ()2.45p in the form 

roc 

M 2)n (zi,ui]z 2 ,u 2 ) = \ dvK(ui,v;zi)M^l(zi,v;z 2 ,u 2 ) + 

JO 

V( Zl , Ul ; z 2 ,u 2 ) + r^(z 1 ,u 1 ;z 2 ,u 2 ) + 0(ra~ 1/2 ), (2.46) 
|r^ 4) (2:i,Mi; 2; 2 ,u 2 )| < q(ui, u 2 )e U2 n~ 1/2 

with polynomial g. The inequality 

\K(u, v-z)\< u l ' 2 v- l l 2 e^ zv (2.47) 

implies that there exists Mq > 2 such that for all z with ?ftz > Mq the norm of the integral 
operator K in the Banach space TL (see (jl.lOj) ) satisfy the inequality 

11*11 < \ (2-48) 

and so there exists the inverse operator (/ — X) -1 . But the problem is that the bound for 
ri above does not allow us to conclude that rff £ TL (recall that we fixe u 2 and consider 
as a function of ui). This difficulty can be easily overcome if we consider a new function 

M 2 , n (z 1 ,u 1 ;z 2 ,u 2 ) = M 2;n (zi,ur,z 2 ,u 2 ) - r ( ^\z 1 , ur, z 2 , u 2 ). 
Then (pT4Uj) takes the form 



M 2>n (z 1 ,u 1 ;z 2 ,u 2 ) = / dvK(u 1 ,v;z 1 )M^l(zi,v;z 2 ,u 2 ) + 

Jo 

V(z 1: «i; z 2 , u 2 ) + K{r^)(z u u x ; za, « 2 ) + 0(^ 1/2 ), (2.49) 
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and (|Z3^) yields 



\K(r^)(z 1 ,u 1] z 2 ,u 2 )\ <C^n- x ' 2 . 



Thus we can apply the (/ — K) 1 to (|2,49p and obtain that for any z : > Mq there exists 
the limit 

roo 

M 2 (zi,ui;z 2 ,u 2 ) := (I — K)~ l (ui,v; zi)V(zi,v; z 2 ,u 2 ) dv. (2.50) 
Jo 

But according to Lemma [T] M 2<n (zi, u\\ z 2 , u 2 ) is an analytic function bounded uniformly in 
each compact in the right half plane of C. Hence, taking any bounded domain U which con- 
tains some z: > Mo, for any fixed ui, u 2 we can choose a subsequence M 2<nk (zi,m; z 2 ,u 2 ) 
which converges uniformly in z\ G U to some analytic in U function. But since for z: > Mq 
for any convergent subsequence there exists a unique limit of M 2>rik (zi, u\ \ z 2 ,u 2 ), defined by 
(|2,50p , on the basis of the uniqueness theorem we conclude that for any z £ U there exists a 
limit of M 2>n (zi, u\\ z 2 , u 2 ) and this limit for $lz > Mq is defined by (|2.50p . Hence we have 
proved (|1.17p for m = 2. 

For arbitrary m we have instead of (|2.4ip 



M n 



n 



/2+l/2 E (o-uiGniz!) 



o(l) 



n n-'D^iz^uj) + f n ( Zj , Uj ) \+0(n- 1 



i=2 



n 



■m/2+l/2 E ( 



o-«iGn(2i) o(l) 
e 11 fn V z h u i) 

3=2 



3=2 



+0{n- 1 ' 2 ) =: Ti + Y, T V + ^(n- 1 / 2 ) 

3=2 

Then, similarly to (|2,43|) . we write T\ from the r.h.s. of (|2.51|) as 
2i = 



(2.51) 



cfo.fir„ (ui, Z\) M myn (zi 
+r%\z\,u\\ ■ ■ ■ ;z m ,u m ) + 0(n~ 1/2 q(ui, . . . ,u m )), 



(2.52) 



where r n admits the bound (|2.44p . 

Since fn does not depend on {aij}™ =2 we can average with respect to these variables 
and, using (|2.37|) write T 2 j in the form 



f 2j 



A o{1) 1 f A o(1) 1 

• 11 fn ( z i,Uj)> =E{V n (z 1 ,u 1 ;z j ,u j )}E< [[ f n (z j ,u j )> 

i=2,i& J l i=2,i^3 J 

+E<V n (z 1 ,u 1 ;z j ,u j ) 11 / n (%,Wj)r 

I „_n i-LA ' 



(2.53) 



16 



Using the Schwartz inequality and Lemmas 1,2, it is easy to obtain that the last term in the 
r.h.s. of (I2T53|) is 0(n -1 / 2 ). Hence, (I23TD . (12331) and (12331) yield 



^m,n \%1 1 ^1 j • • • j -2m; ^m) 



dvK n (ui,v; z 1 )M mtn (z 1 ,v; . . .;z m ,u m ) 



(zi, ui; 2j, Uj)}M m -2 

+0(n- 1 / 2 g 1 ( Ul , . . . ,n m )(e n2 + • • • + e u 0), (2-54) 

Then, using once more the argument, which we applied to (|2.50p . we can prove (|l,17p first 
for $tz > 2 and then extend it to the whole right half plane of C. □ 

Proof of Theorem [7] We prove Theorem U] in two steps: first for polynomial p and then 
extend the statement to any real valued functions p>, satisfying conditions of the theorem. For 
polynomial p we replace in Theorem [3] the product of traces of resolvent of A with different 
Zj (see (I1.19D ) by the product of traces of ipi(A), . . . , p p (A) with tpx, . . . , p m being some fixed 
polynomials. More precisely, we consider (cf (|1.19p ) 



M pn (<Pl, ■ ■ -,Pm) ■= n 



- m / 2 E|Tr p x {A) . . . Tr p° m {A)} = e| fln^Nn^ 



and prove that for any m and any fixed polynomial pi, ... , p m there exists the limit 

lim M m J(pi, ... , p m ) = M m ((pi, ... , p m ) (2.55) 

n— >oo 

and 

m 

M m {p Xl . . . , tp m ) = M 2 (<pi,ipj)M m - 2 (<P2, • • • , <Pj-l, Pj+1, Pm)- (2.56) 

i=2 

Then taking p\ = ■ ■ ■ = p m = P we obtain that there exist the limits of all moments of 

o 

n~ 1 < 2 N n [P] and these moments are expressed in terms of the second moment by the same 
way as for the Gaussian random variable. 

Recall that Theorem [3] imply that the (|2.55p and (|2.56p are valid for tp z (A) = (iX — Zj) . 
We will replace <p z . by the polynomial (fj in (|2.55p - f)2.56|) step by step, starting from the 
last one p Zm (X). To this end we prove by induction with respect to the polynomial degree 
k that if we replace p Zm (X) by a polynomial Pk(X) of degree not exceeding k, then (12.55P - 
(|2.56p are valid. 

o 

For k = 0, 1 N n [Pk] = (recall that Ajj = 0), so ()2.55p - (|2.56p are trivial. Let us assume 
that that we know (|2.55p - (I2.56|) for p m {\) = Pi(\) with I < k — 1 and prove that they are 
valid for I = k. Consider 



Pm{X) = p(X,z m ,k) = -z m X k (i\ - z m ) 1 

= -{-i) k z m (z k m (iX - z m y x + Y, C l k (iX 
\ i=i 



(2.57) 



_ \i-l„k-l 



By the above representation and the induction assumption (|2.55p and (|2.56p are valid for 
Pm(X) = p(X-, z m, k) with any z m . Moreover, if we use the inequalities 



E 
E 



n-^M n [P* k ] 
n- 1/2 N n M\ 



< C(m,k), P fe *(A) = A fc , (2.58) 

< C(m,k)/\%z\, <p* k (X) = X k (iX- z)~\ /c, m £ N, 
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combined with the Holder inequality 

m s 

|M m , n (^, . . .,<p m )\ < [JE 1 /" 1 In-^AfnfcN 

j=l L 

then, since -P^(A) — ip(\; z m , k) = — i^ +1 (A), we obtain 

C 

\M m , n (pi, ■ ■ ■ ,-Pjfc) -M mi „((^i,...,v3(.;z m ,/c))| = | M m ,„ . . . ,(ft +1 )\ < — \ijm ' ( 2 - 59 ) 

I JtZm | 

where C does not depend on n and z m . We will prove (|2.58|) later. Now let us use a simple 
proposition 

Proposition 3 Let the sequence of the functions {n n (^)}^ 1 converges point-wise to the 
function u(z), as n — » oo, in the domain > C , and for any fixed n u n (Q — > u* ; as 

K(C) - <l < Cb/I^Cr, «>0. (2.60) 



T/ten £/iere exist £/te limits 



lim < = lim u(C) (2.61) 



Proof. Take any e > and £ £ such that Co/\$tCe\ a < e/4. Moreover, choose N such that 
\ u n(Ce) ~ u (Ce)\ < £ /4 for any n > N. Then for any n,n' > N 

|< - <'| < K - «n(Ce)| + K' - u n'((e)\ + Wn(Ce) ~ u(Ce)\ + |<v(Ce) ~ ^(Ce)l < £■ 

Hence, there exists it* = lim n _ >0O it* . In addition, for any £ and any e > one can choose iV 
such that |«iv(C) ~~ U (C)\ < e /2 arid \u* N — u*\ < e/2. Then 

NO - u*\ < |«(c) - «iv(C)l + MO - <rl + K + c /l»cr- 

Thus, there exists the second limit in (|2.6ip and it coincides with 

Now if for fixed Zi, . . . , z m _i we consider the functions u n (z m ) = M m , n {ipi, . . . , (p(. ; z m , A;)), 
then (|2.57p gives the point-wise convergence of u n (z m ) and (|2.59p coincides with (|2.60p of 
Proposition [3] with u* = Mk, n (<Pli ■ ■ ■ > P% )■ Applying the proposition we obtain that (12.55P - 
()2.56p are valid if we replace the last function ip m by any polynomial of degree k. 

Repeating the above procedure we replace step by step all (pi,... fm—l by polynomials 
of any fixed degree. As it was mentioned about this implies that for any polynomial P 

o 

n~ 1 ' 2 J\f n [P] converges in distribution to a gaussian random variable with zero mean and the 
variance from (|2.62p . Hence, by the standard argument we conclude that uniformly in x 
varying in any compact of R 

E L™n-VW n [P]\ = e -xy2V(P)^ yr p) = Um VarlrT^AUP]}. (2.62) 

I I n— >oo 

To finish the proof of CLT for polynomials we are left to prove (|2.58p . It is done in the further 
proof of Theorem [H 

To extend CLT to a wider class of functions we use 

n 

Proposition 4 Let be a triangular array of random variables, M n [(p] = vi^i^) 

l=i 

be its linear statistics, corresponding to a test function ip : R — > R, and 

KM = Var{n~ 1/2 .A/" n [<£>]} 
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be the variance of M n [p\. Assume that 

(a) there exists a vector space C endowed with a norm ||...|| and such that V n is defined 
on C and admits the bound 

V n [ip\ < C||v;| | 2 , yip G C, (2.63) 

where C does not depend on n; 

(b) there exists a dense linear manifold L\ C C such that the Central Limit Theorem is 
valid for M n [p], p> G £\, i.e., if 



Z n [xtp] = E < e 



Jxn 1 / 2 Af„[ip] 



is the characteristic function of n~ 1 ' 2 Af n [ip], then there exists a continuous quadratic func- 
tional V : C\ — > K+ such that we have uniformly in x, varying on any compact interval 



lim Z n [x(p] = e 



V</? G C\\ 



(2.64) 



Then V admits a continuous extension to C and Central Limit Theorem is valid for all N n [<p] , 
<pe£. 

Proof. Let {pk} be a sequence of elements of C\ converging to cp G C We have then in 

o 

view of the inequality \e m — e lh \ < \a — b\, the linearity of M n [(f] in ip, the Schwarz inequality, 
and (12331) : 



Z n (x(p) - Z n (xip)\ 



V =l Pk 



< HE 



n-^Mnl^-n-^Mnipk] 



< IxlVsLr^in-^Nnitp - <p k }} < C\ 



x\ \\P - <Pk\ 



Now, passing first to the limit n — > oo and then k —* oo, we obtain the assertion. □ 

Let us show now that hypothesis (a) and (b) of Proposition U] are fulfilled in some vector 
space. We fix some c > and consider the vector space C of functions <p such that (p, cp', ip" G 
L 2 (M,cosh- 2 (cA)) (see (fl~22]) ). Denote 



p"{\)\ 2 cosh' 2 {c\)d\+ / \p'{\)\ 2 cosh' 2 (c\)d\+ / \^(X)\ 2 cosh- 2 {c\)d\ 



It is evident that the space of all polynomials is dense subspace in £ with respect to 
the norm ||.||. Moreover, (|2.62p proves (b). Hence we are left to check assumption (a) of 
Proposition [H 

It is easy to see that if p G C then /(A) = p>(\) cosh _1 (cA) G L 2 (M) and also /', /" G L 2 (M) 

and 

ii/iiL w + iiriiL(R)<cibii 2 . 

Hence it is enough to check that 

Vzrln-^Tr f(A)e ±cA } <C(\\f\\] 
According to Proposition [T] (see (|2.4|) ) 



e " 1 1 2 



(2.65) 



E 



■WTr f{A)e ±cA 



2m 



< C2 m E 



Tr (f(A)e ±cA -f{A^)e 



2 m 



(2.66) 
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where ^l^ 1 -* is defined in (j2.6f) . Note that to prove (|2.65p it suffices to consider m = 1, but we 
need other m to prove (|2.58p . Write 



Tr (f(A)e 



cA 



f(AW)e 



where / is the Fourier transform of /. Then the Duhamel formula yields 



Tr ( e (i?+c)A - e (i?+c)A(1) 



£ dt{it + c)Tr (e^+ c ) A {A - A^)e^ AW ^ 

/ ^ + c) f>^ +c)A ),i 
Jo ~\ 



f e (ie+c)A(i)(i_ t ) a (i)> 



2tcA 



ei,ei 



| 1 /2 e 2cAW(l- t ) a (l) )a (l) ) l/2 



Here we used that 

Tr (e^)A {A _ A( i ))e(i z +c)A W(i-t^ = ^( e *«+c)A )iiaifc(e «+c)A(«(i-t) )w 

The first some gives ^™ =1 (e'^ +c ) A ei) i (e^ +c ^ (1) ( 1 -*)a( 1 )) i where e x = (1,0,... ,0) and the 
vector o« is defined in (|2.8jl . The second sum is 0, since relations = A$ = (i = 
1, . . . , n) imply (e^ +c ^^^ 1— = (i = 1, . . . , re). Then the Schwarz inequality yields 

Var{re- x / 2 Tr f{A)e cA } < C 2 (J + c^Ve 1 /^^, ei)2} 

EV2{( e 2(i-^AM a (i)yi) )2} _ (2>67) 

Using the Schwarz inequality once more and then the symmetry of the problem, we obtain 

E{(e 2tcA ei,ei) 2 } < E{(e itcA ei , ei)} = E{Tr e 4teA }. 

Similarly, using the Schwarz inequality and then the independence .A^ 1 ) of aft) , we can average 
with respect to to obtain 

< C^+p^EK^r e 4(l-*)cAW } + c(p 2 +p 3 )E l/2 {n -l Tr e 8(l-t)MW } _ 

Since all entries of A and A^ and A — A^ are positive, we have for any t 

E{Tr e 4(i-*M (1) } < E{Tr e^M} < E{Tr e 4cA }. 
Moreover, according to the result of [5] we have for any m 

E{n- lr Tr A 2m } < C^ml =► E{n _1 T> e cA } < 2e c ° c2/2 . 
In addition the Schwarz inequality yields 

l/(0l(iei + <W) < / i/(0i 2 (iei + c) 4 ^ /(iei + c)- 2 de<c(n/n 



lL 2 ( 



+ 



f " 1 1 2 

IIl 2 ( 



Summarizing the above inequalities, we obtain (12.650 and hence the assumption (a) of Propo- 
sition UJ Then Theorem U] follows from Proposition 01 

To prove (|2.58|) we use again (I2T66D . where for the first line of QTMb /(A) = A fc cosh -1 (cA) 
and for the second line /(A) = X k cosh _1 (cA). Repeating the above argument we obtain 
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